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Abstract 

The unique real root of cos(a;) = x, recently referred to as the Dottie 
number, is expressed as an iteral of cosine. Using the derivatives of iterals, 
it is shown why this number is achieved starting from any real number, 
when the iterates of cosine successfully approach infinity, and how this 
affects the Maclaurin series of the iterals. Properties of the iterals of 
cosine and sine and their derivatives are considered. A CH — h template for 
iteral is applied for computation of Julia sets. 

1 Introduction 

If one sets a calculator to work with radians, types any number, and presses the 
button cos over and over again, then soon the screen shows unchanging further 
decimal digits 0.739085133215.... This is a truncated root of the equation 
cos(a;) = X. Remembering such a story occurred with Dottie, a professor of 
French, Samuel Kaplan refers to the root as the Dottie number [7|. He formulates 
exercises involving the number for Advanced Calculus, Differential Equations, 
Chaos, Complex Variables, and concludes 

It is unlikely that the Dottie number will enter the annals of great constants 
along side e, it, the Golden Ratio and many others. However, the Dottie number 
and its story might make good teaching elements for others out there. I also 
imagine there are many other interesting facets of the Dottie number yet to be 
discovered. 

The two properties that this is 1) the only real root of cos(a;) — x, and 2) 
the attr actor in the iterates of cos, where the domain of attraction is the entire 
real line had been known prior to the Kaplan's naming proposal. 

The first property implies computation of the root. In 19th century, T. Hugh 
Miller |S] evaluates eight decimal places by a method of trial. He mentions that 
the problems of finding the angles satisfying the equation occur in analytical 
geometry and presents a proof of the uniqueness of the real root 

For while x increases from to ^, cos(a;) diminishes from 1 to 0. Also cos(a;) 
is positive while x lies between and —1. 

The chart of the line y = x and curve y — cos(a;) intersecting only in one 
point is a geometrical proof. Then, for this root cos(a;) — x = arccos(a;). The 
first hundred decimal digits of the Dottie number are found in [7] and more in 
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[5]. Hrant Arakelian presents 6,400,000 of them [T]. Miller proposes a method 
for computation of complex roots ^ and reminds that their number is infinite 
[To] . In [9] he notices without a proof 

It is easy to show that one value of the angle is cos cos cos . . . cos 1, when the 
cosine is taken successively to infinity. 

He does not discuss iterative properties of cosine. It is not clear, whether 
Miller knew that "cosine taken successfully to infinity" gives the same result 
for any real number. The Miller number is a suitable name, if the community 
wants to emphasize his computational contribution. 

The second property is discussed, for instance, by Robert Devaney jlj pp. 
39 - 56]. Applying graphical analysis, he illustrates dynamical properties of the 
functions ^/x, a;^, 4x(l — x),x^ ~2, sin(a;), and cos(x). Graphical analysis shows 
that iterates of cosine starting from arbitrary real values tend to the real root of 
cos(a;) = X. This method is used by researchers to study function iterates. For 
instance, Mitchel Feigenbaum discovers a number named today after him [S]. 

Probably, generations using calculators have faced with the "remarkable cos 
experiment". Psychology of "unexpected" is responsible for the first impression. 
Many would agree without a calculator, that successive iterations of y/x tends 
to one beginning from any x > I and, after more thinking, that this is true for 
X > 0. But for cos it is not obvious. One is easier to memorize than 0.739085 . . .. 

Since 1980th, Hrant Arakelian emphasizes the fundamental meaning of the 
real root of cos(x) = x viewing it as a hidden link of mathematics and suggest- 
ing the names cosine constant, cosine superposition constant, cosine attractor, 
cosine fixed point [T] . He proposes to denote it by the first lower case character 
of the Armenian alphabet pronounced [a:] and claims relationships between this 
and some fundamental physical constants. If his findings will be eventually re- 
cognized, then the name Arakelian number will become inevitable. The author 
is not yet ready to share Arakelian's views. 

It is natural for the author jTT] to apply to the iterates of cos(a;) the notation 
for the iteration of functions recently introduced and referred by him to as the 
iteral 

Mj, cos(a:;). 

For a function of one variable, v is the initial value of the variable, and n is the 
number of the iterations. On each iteration a previous function value becomes 
its new argument. By definition, the iterals of cos(x) are 

cos(a;) — V, M„ cos(x) = cos(w), M„ cos(a;) — cos(cos('y)). 
The "intriguing calculator experiment" hints that 

vw e R Mr cos(a;) — Dottie number = D. 

Here, V "for all" is the universal quantification of proposition, and R denotes 
the set of real numbers. The goal of this article is to review the properties of 
the iterals of cosine and sine, derivatives of the iterals, Maclaurin series, and 
show why the Dottie number is achieved. 
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2 Iterals of Cosine and Sine and First Derivatives 



By the definitions of the Dottie number and iteral 

Vn € No M cos(x) — D. 

No is the set of natural numbers and zero. This reflects that Z? is a fixed point 
of cosine. Since cosine is even function, the following bounds apply 

Vw e R M y cos(x) e [-1, 1] and cos(a;) e [cos(l), 1]. 
If the function values belong to its domain up to n iterations, then 

For cosine and sine this is the case for all u € C and n > 0, where C denotes the 
set of complex numbers. Iterals of cosine and sine are the functions of the initial 
value V. They create two families of functions, where each member corresponds 
to the order of iteral n. When the initial value w is a valid expression containing 
the variable x, parameters, and functions the iteral creates new functions of x. 
The function under the iteral sign can be a vector function 

Replacing v with x makes the iterals of cosine and sine functions of x 

cos(x) and sin(x), 

For basic theorems on iterated functions see [2] . Let us find their first derivatives 
using the chain rule and iteral notation. We get for the orders n = 0, 1, 2 

n cos(x) (J^x cos(x)^ sin(x) (^^x sin(2;)^ 

Ox 1 X 1 

1 cos(a;) — sin(a;) sin(a;) cos(a;) 

2 cos(cos(a;)) sin(a;) sin(cos(a;)) sin(sin(x)) cos(a;) cos(sin(2;)) 



With the increasing order of the iteral of cosine the previous first derivative 
multiplied by — sin(M^ cos(a;)) yields the current first derivative. For the 
first derivative of the iteral of sine the multiplier is cos(Mj. sin(x)). Thus, 



= (-1)" n sin(M!cos(x)) and 

fe=0 

dM^sin(x) T/r'= . 
" =11 cos(yi,sm(x)). 



dx 

k=Q 



For rt = no multiplications are performed. The default value of Y[ is one. The 
iterals make the expressions compact for large n. Let us prove the formula for 
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the iterals of cosine using mathematical induction. The proof for the iterals of 
sine is similar and omitted. The formula is correct for the orders 0, 1, 2. Let it 
be valid for n > 2. Then, for rt + 1 we get 



dH"^' COS(X) _ ^^HIcos(x) ^°^(^) _ dcOs(H"cOs(x)) _ 

dx dx dx 

= -sm(M:cos(x))^^%^ = 
= - sin(M, cos(2;))(-l)" J]^ sin(M, cos(2;)) = [| sin(M, cos(a;)) = 

k=0 k=0 
h-1 J, 

= {-1)'^1[MV[^cos{x)). 

k=0 

This is exactly the same formula for n replaced with h = n + 1. 

Since for x G [—1, 1] sine is a monotonously increasing function and — 1 < 
-0.84... = sin(-l) < sin(a;) < sin(l) = 0.84... < 1, each factor, with the 
exception of the first one, in the first derivative of a cosine iteral of the order 
n > 1 is by absolute value less than 0.85. This is because of the mentioned 
bounds of the cosine iterals serving as arguments for the sines. For any x the 
iterals of the order n form a variant, where each member starting from the 
second by absolute value is less than \q = 0.9|". The latter is infinitesimal 
variant for n — >■ oo [6, p. 50]. Then, 



dK 



Vx e R lim "f"''^'^^ = lim (-1)" ff sin(M, cos(a;)) = 0. 

fe=0 

Because all iterals of cosine also intersect at the fixed point D this ensures that 

Wx e R lim cos(a::) = D. 

n— voo 

The idea to consider the limit of the first bound derivative of the iterals suggested 
by Alexander Tumanov in a private discussion simplifies the proof. The latter 
formulated by the author with a variant of iterals complements the Graphical 
Analysis proof and Miller's notice about cos cos cos . . . cos 1 iterated to infinity. 

For n > iterals of cosine are even functions, below N is the set of natural 
numbers without zero, 

_^ cos(x) = cos(a;). 
The iteral of the order zero x is odd function. All iterals of sine are odd functions 
Vn e No M_3, sin(a:) = -M^ sin(a;). 
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For n > the iterals of cosine and sine are periodic functions. This follows 
from a theorem on iterates of periodic functions [2]. For n = 1 we deals with 
cos(x) and sin(a:) and their period 27r. For n > 2 the period of iterals of 
cosine is tt. This is because cos(cos(a; + vr)) = cos(— cos(a::)) = cos(cos(2;)). 
Further composition by cos cannot change the result obtained after the two 
initial iterations. Increasing the order of iterals of sine does not reduce the 
period for n > 2 because sin is odd function. It remains equal to 2n. 

Independently on n, the real domain of each iteral of cosine and sine is 
(— oo, oo). The ranges of the iterals for n = are also (— oo, oo). The ranges of 
the iterals for n = \ are [—1, 1]. Again, this is true for real functions of cosine 
and sine. For n = 2 the ranges of the iteral of cosine and sine are [cos(l), 1] and 
[— sin(l), sin(l)]. Till now, the ranges are narrowing with the increasing order 
n for both families. This remains true for n > 2. Indeed, on the interval [—1, 1] 
sin is monotonously increasing function. It takes its minimum value at —1 and 
maximum value at 1. The minimum is greater than —1 and the maximum is 
less than 1. They form a narrower input range before the third iteration of sin. 
Since x < sin(a;) for x G [—1,0) and sin(a;) < a; for a; G (0, 1] a current iteral 
range is narrower than a previous one. Since sin, being odd, is antisymmetric 
with respect to the origin of the coordinates, the bounds are moving towards 
zero so that it remains the middle of the iteral range. Thus, Va; e R V n e N 

M!i sin(x) < M!isin(x) < M"sin(x) < M" sin(a;) < sin(x). 

The n order iterals of sine of 1 (-1) form a decreasing (increasing) variant with 
positive pn (negative (/„) members. Therefore, the variant p„ (g„) has a limit: 
non-negative a (non-positive b). Let us consider sin(p„) (sin(g„)). On the one 
hand, its value at the limit is sin(a) (sin(6)). On the other hand, it is also an 
interal of the order n + \. Thus, sin(a) = a (sin(6) — h). This is possible only, if 
a — h = 0. The final part of the proof with "extra" sine applied to the variant 
is given by Alexander Tumanov in a private discussion. Hence, 

Vx e R lim sin(x) = 0. 

For the iterals of cosine the picture is richer. The range of the second order 
iteral is [cos(l), 1]. This becomes the input interval for the last cos evaluated for 
the third order iteral. On this interval cos is monotonously decreasing and posit- 
ive function. Accordingly, it takes its minimum cos(l) and maximum cos(cos(l)) 
values at 1 and cos(l). The range of the third order iteral is [cos(l), cos(cos(l))]. 
Let us notice that this range is narrower than for the second order iteral because 
cos(cos(l)) < 1 and the low bounds for the second and third orders are identical. 
By these reasons, cos on this interval is also monotonously decreasing and posit- 
ive. Meanwhile, this range becomes the input interval for the last cos evaluated 
for the fourth order iteral. It takes its minimum cos(cos(cos(l))) at the most 
right argument cos(cos(l)) and its maximum coscos(l)) at the left most argu- 
ment cos(l). The range of the fourth order iteral is [cos(cos(cos(l))), cos cos(l))]. 
Let us notice that this range is narrower than for the third order iteral because 
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cos(l) < cos(cos(cos(l))) and, now, the top bounds for the third and fourth or- 
ders are identical. Then, on this input interval for the last cos of the fifth order 
iteral cos remains monotonously decreasing and positive. During two consecut- 
ive steps cither low or bottom bounds arc identical and then move towards each 
other. On any closed subinterval of [cos(l), 1] cos evaluated at the most right 
end is less than cos evaluated at the most left end. A range of the order n > 1 
iteral of cosine can be expressed as 

rn-l-(n%2) jrn-2+{n%2) 



I cos(a;), cos(x) 



where %2 is division modulo 2. The three sequences clarify the formulas 

n23456789... 

n-1- (n%2) 1 1 3 3 5 5 7 7... 
n-2 + (n%2) 02244668... 

yielding 

Hn-l-(n%2) Ti-Ti -r i-Ti-2+(n%2) 

1 cos(a;) < Fi^cos(a;) < Fii cos(a;). 

Being iterals of cosine, the bounds tend to the Dottie number with n — >■ oo. 

The locations of minimums and maximums of the iterals of cosine are the 
roots of ^ 

cosfa;) , tt • /T/T*^ / x\ 
^^^^ = (-1) n cos(ar)) = 0. 

fe=0 

From all sin factors in the product only sin(a;) and sin(cos(x)) can get zero 
values. For n > 1 the iterals of cosine in the remaining sin factors have the ranges 
within the interval [cos(l), 1]. On this interval sin(a;) > 0. Thus, the minimums 
and maximums of the n > 1 iterals of cosine at the roots of sin(.T) sin(cos(3;)) = 
are equal to ±!^, k gNq. 

The locations of minimums and max;imums of the iterals of sine are the roots 

of 

!^%W = ncos(H!s.„W)=0. 

fe=0 

For the iterals of sine of the order n > each expression gets the factor cos(a;). 
Therefore, the locations coincide with those of sine because these iterals have 
the same period as sine. They are equal to | ± kw, k G Nq 

The following figures depict the first seven iterals of cosine and sine, locations 
and sizes of maximums and minimums and intersection points. The unbound 
iteral of the order zero is the straight line y = x for both cosine and sine. It 
is omitted in order to maintain suitable drawing scale. Due to evenness and 
periodicity of the iterals of cosine, all intersection points correspond to y = D. 
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D is Dottie number 1 cosfx) 

2 cosfcosfx)) 

3 cosfcosrcosM)) 

4 cosfcosfcosrcosfx)))) 

5 cosfcos(cos(cosicos(x))))) 

B cosfcosrcosfcosrcosrcosm))))) 

7 cos(cos(cos(cos(cos(cos(cos(x))))))) 




1 sin(x) 

2 sinfsinfx)) 

3 sinisin(sm(x))) 

4 sinisinfsinrsinfx)))) 

5 sinisinisinfsinfsinrx))))) 

B sinfsinfsinfsinfsinfsmfx)))))) 

7 sin(sin(sin(sin(sin(sin(sin(x))))))) 
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The alternating distances between the intersection points on the cosine curve, 
iteral of the first order, are h — 2D and Li = 2(7r — D). The alternating 
distances between the intersection points on the higher order n > 1 iteral curves 
are ^„ = 2D and L„ = tt — 2D, where | > 2D. Here, li = 1^ and Li — 2Ln + ln. 
The line y = D is never the center of the iteral of cosine ranges. This is the 
cause of the alternation. 



3 The Series of Iterated Cosine and Sine 

Using D'Alembert's ratio test it is easy to establish that both cosine and sine 
Maclaurin series 

^3 5 °° ^2n+l 

sin(z) = z-- + --... = V(-ir 



3! 5! ' (2n + l)! 

n— 



are absolutely convergent for any complex argument z. Thus, the series repres- 
enting the second power of cosine and sine expressed as the Cauchy products 
also converge absolutely 



2n \ , f n \ n\ 

I 



cos\z)^y{~lTi—y[ , where ^ =77^^ 



CX3 



sin^(z)^y(-ir ^ 



' (2n + l)! V 2fc + l 

n=0 ^ ^ fc=0 ^ 

Here, cos^(2) and sin^(z) are the ordinary second powers of cosine and sine. 
Due to the new iteral notation they cannot be misinterpreted in mixed 
expressions as the twice iterated trigonometric functions. 

The Cauchy products sequentially created for higher powers of each function 
from the absolutely convergent components are absolutely convergent series on 
each step [S]. Our interest is even powers of cosine and odd powers of sine. 

One way to get the series for cos(z) = cos(cos(z)) is to replace z in each 
power term of the Maclaurin series expressing cosine with the series itself. For 
instance, the first two terms are 

M; cos(z) = cos(cos(z)) = Yll Y.{-ir^ ^ 1 - + . . . = 



n=Q\ ^) (2n)! ^k=0 y 2k J 



2! ' "^"^"^2! " 2[^32['"-'^ 
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Due to the Cauchy product construction, for any power of cosine there is one 
term in the series not containing z and equal to one. Combining these units 
divided by the corresponding factorials gives for the second order iteral of cosine 
the infinite scries l^^ + 2T^gT + ^^''' followed by the terms containing 
the even powers of z. Absolute convergence ensures that regrouping the terms 
retains the limit value. The initial series independent on z is equal to cos(l). 

All powers of z coming from the original series for cosine will stay in the 
Cauchy products. Presence of units prevents their vanishing. For instance, 
the second powers of z found in the increasing powers of cos^"(2;) come for 
n = 1, 2, 3... as — z^, — 2z^, — 3z^, ... — nz^. After replacing z in the series for 
cos with the same series, each of these terms is multiplied by its own factor 
(-1)" 2^ yielding 

cos{z) = cos(cos(2;)) = cos(l) + - ^ + ^ )z'^ + -- - ■ 

Multiplying and dividing the terms within the intermediate brackets by two 
gives i(|-| + |-... + (-l)"+i^-...) = isin(l)and 

Mj, cos(2;) = cos(cos(2;)) = cos(l) + ^ sin(l)2:^ + higher even powers of z. 

Let us compare this series with the series for the iteral of cosine of the order one 

Mj, cos(.2) = 003(2;) = 1 — ^z^ + higher even powers of z. 

We notice that a) the first term independent on z in the iteral of the second order 
is the first iterate of cos of the first term of the iteral of the first order; b) the 
new coefficient at the power z'^ is less by absolute value than the corresponding 
coefficient in the iteral of the first order; c) the sign of this coefficient alternates. 
The successive series may not get odd powers of z. While with the increasing 
number of iterations the coefficients containing powers of z are less obvious at 
the moment, the term not containing z is getting clear 

Mj, cos(2:) = cos{z) + terms containing z. 

Another way to the same conclusion is representing an iteral of cosine of the 
order n by the Maclaurin series (compare it with the previous series) 
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H:cos(.) ^ M:cos(.) + ^^%^(.o - 0)|, + '^%^(.o - 0)|, 

= nr cos(.) + (-1)" n -(m: cos(.)) ^ + ' - 



(^0 = 0)- + ... 

fe=0 



where the first term is f{zo = 0) = Mg cos(2;) = cos(2;). Reduction of 

the order by one takes place because of the first iteration cos(2;o = 0) = 1. Based 
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on the Miller's notice and our (Tunianov and author) proof 



lim Mj^ cos(2;) = Mj^ cos(z) = D. 



Since the limit of the iterals of cosine as it was proved does not depend on the 
initial real value, the limit of the remaining terms exists and is equal to zero 



cos(z) = cos(2:) + lim terms containing z = D 



0. 



The "calculator experiment" confirms it for the real initial values. Computer 
experiments show that for some initial complex z the iterates of cosine approach 
the Dottie number. However, for some other the real and imaginary parts grow 
and cause the overflow condition. While the expression for the first derivative 
of the iterals of cosine and sine as well as the series constructed in this section 
remain valid for the complex arguments, the limit on the iteral order does not 
have to hold because successive iterals and the first derivatives are unbound for 
some complex arguments. To clarify details we need expressions for higher order 
derivatives of the iterals of cosine. A simplifying condition is that all derivatives 
for Maclaurin series are evaluated at zq = 0- 

4 Higher Order Derivatives 

The first derivative of the iteral of cos of the order m > can be written as 

jT/T"^ / V m — 1 ^ m— 1 ^. 

dz " n sin(M,cos(z)) (-l)"sin(z) [] sin(M, cos(z)) 

fc=0 fe=l 

It is equal to zero for z = because of the factor sin(2 = 0). Thus, the term 
/ (zo-o) ^ vanishes. In the limit case m — cx) we reuse the property that the 
infinite product is equal to zero, if at least one of its factors is equal to zero. 

The second derivative is taken as the first derivative of the product of sines 
composing cosine iterals of increasing orders. In courses of Calculus, it is proved 
using mathematical induction that [uvw...s\' = u'vw...s + uv'w...s + uvw' ...s + 
uvw.-.s' in p. 201]. The author writes the formulas involving the number of 
factoring functions fi{x),i — 1, ...,m. If fl{x) exists for all i, then 

(m \' rn ( / i-1 \ / ™ ^ 

i=l / i=l [ y=l / 

If we put jrj^ — 1 also when fi{x) = and take care grouping and computing 
this ratio first, then simpler formulas are 
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Each formula in the right side "skips" the factor fi{x) in the ith summand. We 
notice that in the last formula = ln(/i(x))'. Applying the last substitution 
implies that fi{x) ^ 0. All summands but the first one have the factor sin(z). 
They are equal to zero at Zq = 0. The non-vanishing term is 

' ^ ' {zo = 0) = (-l)'"cos(^o = 0) n sm(Fio cos(^)) = 

^ fe=i 

= (-1)™ n sm(l/lo cos{z)) = (-1)"' Y[ sin(l/li cos(^)). 

k=l fe=l 

For m = 2 we get sin(l). This completes the coefficient for z"^ in the Maclaurin 

series of the second order itcral of cosine as | sin(l) previously obtained by a 
different method. Let us notice that 

{-irl[sm{ll, cos(^))j|j =0. 

This is because the iteral of cosine with the initial value one tends to D with 
m — >■ oo (Miller's notice). Due to the inequalities 0.6 < sin(£)) < D < 0.8, the 

2 

product forms an infinitesimal variant multiplied by ^ selected in advance and 
playing the role of a constant. 

In order to move forward we need a formula for nth derivative of the product 

of m functions, where each function has nth derivative. For m > and n = 1 
the formulas are constructed. For m — 2 and n > this is the Leibniz formula 

(u(a;)w(a;))(") = ( " ) u^''-'\x)v^'\x). 

Here, the "powers" in brackets denote derivatives of the corresponding orders. 
The order zero means the function itself. Using binomial, trinomial, n-nomial 
trees corresponding to the second, third, nth derivative it is easy to show 
that the sum representing the nth derivative of the product of m functions will 
contain m" summands. Each is a product of m factors. Each factor is either 
a derivative of a certain order or a function. If the function is viewed as a 
derivative of the order zero, then the sum of derivative orders for each product 
must be equal to n. Some of the products will be identical and can be combined 
getting the mTiltinomial contributing coefficients. This will reduce the number 
of summands. However, the sum of the coefficients still must be equal to m". 
Then, for i = 1, 2, . . . , m and ji = 0,1, ... ,n 

(m \ (") / , m \ 



While each index ji can get any value from to n, summation is done on all 
those combinations of the indexes, where the sum of their values remains equal 
to n. Let us verify that the formula works in a few already known to us cases. 
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If m = 2 and n = 0, then it reduces to the product of two functions. In this 
case the only summation combination is (ji = 0, j2 = 0) and the only coefficient 
is equal to ^ = 1. Then, (/i(a;)/2(a:))^°) = fi{x)f2{x) from the left side is 
equal to lf^\x)f^\x) — fi{x)f2{x) from the right side. 

li m — 2 and n > 0, then j2 = n — ji and the relevant n + 1 combinations are 
{ji ,n — ji). The coefficients are equal to ^. jj^ , = ^-j^^lj^-), = binomial coeffi- 
cient (n,ji). The product of functions contains only two factors f[''''\x)f2"~''^\ 
Thus, the right side reduces to the Leibniz formula. 

If m > 2 and n = 1, then the relevant combinations are those, where one 
index is equal to one and all other are equal to zero. The total number of 
combinations in this case is equal to m. Thus, all coefficients are equal to one. 
Since we included the derivative of the order zero into considerations, the right 
side is the sum of YihLi //"'''' (2;), where only one factor is the first derivative. 
This is exactly our result obtained earlier. The case ni — 1 and n > is 
also covered by the formula. The formula exposes similarity with expansion of 
{xi + ... +Xm)^ as it was pointed to (without giving the formula) in [6, p. 238]. 
Here, n without brackets is the ordinary power. The author views the formula 
as well known and omits the double mathematical induction proof on m and n. 

This formula should be applied to the first derivative of cosine iterals. Due to 
the presence of the factor sin(z) in the first derivative, each summand in the just 
proved formula will get the factor sin(2:o = 0) = for each coefficient in front of 
an odd power i of Thus, all those terms vanish from the Maclaurin series of 
the iterals of cosine. This is understood, because the "last" cos series contains 
only even powers of z and all replacements also contain only even powers of z. 

2p 

The coefficient of the Maclaurin series at for p > 1 is 

m-l ^^ \ (2P-1) / „ N (2p-l) 

nsin(M,cos(z))) ={-ir[l[sm{V[, cos(z)) 

\k=l / 

The power is odd 2p ~ I because we start from the first derivative given as a 
product of functions. The formula for the nth derivative of the product gives 




4)™ nsin(M, cos(z)) 



= (-ir E (l|^n(--(Ml"cos(.)) 

Let us prove that this coefficient tends to zero in the limit m — > 00 for all p > 1 
at z = 0. For p = 1 we have proved it already. 

The first function under the product in the left side is sin(z). Its odd derivat- 
ive 2p—l is equal to cos(z) with a precision of sign. Thus, the summand for the 
combination (ji ~ 2p — 1, 72 = = •■• = jm = 0) remains with the first factor 
cos(2;o = 0) = 1 and the multinomial coefficient equal to one independently on 
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TO. All other factors for I > 1 are sines of iterals of cosine sin(J/lQ cos(z)). 
This construction as we have seen tends to zero with to — >■ oo as an infinite 
product of the numbers closer and closer approaching 0.6 < sin(D) < 0.8. 

If for / > 1 at least one index ji = 1. then the sine of the iteral of cosine 
being a compound function produces the factor sin(z), which will make the entire 
summand equal to zero at zg = 0. Summands with ji — even number will be 
equal to zero because of the derivative of the first factor sin(z). Summands with 
ji = will be equal to zero because they will contain the factor sin(z) as the 
first function. We see that the question about equality of a summand to zero 
at zq = translates to the question in the number theory and consideration of 
individual > within the constraining sum + + = 2p—l. Here, 

the right side is fixed for a given p, while the left sum gets increasing number of 
non-negative terms for to oo. The infinite "majority" of them must be equal 
to zero leaving the function-factors "untouched". The product of the untouched 
factors will tend to zero at zq = with to — cx). Without losing generality 
(since other cases are already considered) , we set ji — 2q— l,l<q<p and get 
j2 + ---+ji + ---+jm — 2{p—q). The number of "touched" (differentiated) factors is 
finite for any given p. Each of them is bound (they are computed at real zq = 0). 
Thus, the product of a finite number of bound factors and infinite number of 
factors tending to zero, should tend to zero. Presence of multinomial coefficients 
does not infiuence on this conclusion. Thus, applying the Maclaurin series we get 
for the terms containing z in the limit of the iteral order approaching infinity the 
sum of the infinite number of infinitesimal variants. We know that it vanishes 
for real z due to our and other proofs such as graphical analysis or Newton's 
method. For complex z it is not so. 



5 Filled in Julia Sets of Cosine and Sine 



Using the standard set building notation and iterals, a filled in Julia set can be 
defined as z S C, Jc = {v & C : |HJf'/(z)| < oo} [U. The definition of the 
Mandelbrot set is another example of usefulness of the iterals: z g C, M = {c G 
C : |Hg°(z^ -I- c)| < oo} [TT]. For cosine and sine the filled in Julia sets can be 
defined as ^ 

z e C, Jcos = {v e G : Mj, cos(z) < oo} 



e C, Jsin = {w e C : M„ sin(z) < oo} 



Computer experiments show that for some complex numbers iterals tend to 
infinity by absolute value. The structure of the sets is sophisticated. The points 
on a complex plane defining the sets form queer patterns. Appendix presents a 
C++ program and figure of a subset. 



Acknowledgments. I would like to thank Alexander Tumanov for reviewing 
the article, preventing penetration of an error into it, and two ideas completing 
two proofs marked in the text. 
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6 Appendix 



The author has written the CH — h program iteral_juHa_set.cpp. The Standard 
CH — h Library has the class complex and a set of overloaded functions using the 
complex numbers [12 . The comment is self-explanatory. 

// Progrsun scans a rectangular area of the complex plame specified 
//by two points xl + i * yl Eund x2 + i * y2. The grid is the number 
// of points along x- and y-sides of the rectangle. The last parsuneter 
// selects iteral of cos or sin. Parauneters should be given in the 
// command line. Exsunple: 

// iteral_julia_set -2.5 -2.5 2.5 2.5 500 cos > f.txt 

// Progreun outputs coordinates of complex points, where an iteral tends 

// to inifinity. The output format is friendly to Gnuplot progrsun. In 

// the excunple the output is redirected to the file f.txt. This file 

// can be used with gnuplot or wgnuplot applying a single command 

// plot "f.txt" with dots 8 

#include <stdexcept> /* Stsindard C++ Library, exceptions */ 
#include <iostrecun> /* Standard C++ Library, I/O stresuns */ 

#include <iomaiiip> /* Stctndard C++ Library, manipulators */ 

#include <string> /* Stctndard C++ Library, strings */ 

#include <complex> /* Stctndard C++ Library, complex numbers */ 

#include <cmath> /* Stctndard C++ Library, cos etnd sin */ 

using Ucunespace std; 

// The Iteral template for a function of one variable. 
// n is the number of iterations >= 
// V is the initial value 

// f is the function such as cos, sin, etc. 

template<class T, class F> 

T iteral (unsigned int n, const T& v, F f) 

{ 

T tmp = v; 

forCunsigned int i = 0; i < n; ++i) tmp = f(tmp); 
return tmp ; 

> 

int main(int argc , char*argv[]) 
{ 

try { 

if (argc < 7) { 

cout << "Usage: " « argv[0] « " xl yl x2 y2 grid cos I sin" 

<< endl; 
return 0; 

} 

const unsigned int n = atoi(argv[5] ) , ITERATIONS = 50; 
if (n < 2) 

throw invalid_argument("Grid (" + string(argv[5] ) 
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+ ") must be >= 2") ; 
if (string(argv[6]) != "sin" && string(argv [6] ) != "cos") 
throw invalid_arguinent("T3rpe sin or cos but not " 
+ string(argv[6])) ; 
complex<double> (*f ) (const complex<double>&) = 
string (argv [6] ) == "cos" 

? (complex<double> (*) (const complex<double>&) )cos 
: (complex<double> (*) (const complex<double>&))sin; 

const complex<double> zl (atof (argv [1] ) , atof (argv [2] ) ) ; 

const complex<double> z2(atof (argv[3] ) , atof (argv [4] )) ; 

const double dzr = (z2 - zl).real() / (n - 1); 

const coinplex<double> dzi(0.0, (z2 - zl).iinag() / (n - 1)); 

complex<double> z(zl); 

for (unsigned int r = 0; r < n; ++r, z += dzr) { 

for(unsigned int i = 0; i < n; ++i, z += dzi) { 
if (normdteral (ITERATIONS, z, f)) < 10) 

cout « setw(25) « setprecision(16) « z.realO 
« " " « setw(25) « setprecision(16) 
« z.imagO « endl; 

} 

z = coinplex<double>(z.real() , zl.imagO); 

} 

} 

catch(const exceptionfe e) { 

cerr « e.whatO « endl; 
return -1; 

} 

catch ( . . . ) { 

cerr « "Unknown exception" « endl; 
return -2; 

} 

return ; 

} 

For each point in a rectangular grid the program evaluates iterals of cosine 
or sine for 50 iterations. If the itcral value is less than 10 by absolute value, 
then corresponding coordinates of the point are output. Otherwise, the point 
is skipped in the output. The standard output can be redirected to a text 
file and plotted using gnuplot or wgnuplot shareware. The result of plotting 
500 * 500 = 250000 scanned and checked points in the square —2.5 — 2.5i and 
2.5 + 2.5i are on figure. Setting ranges and grid of the rectangle works as zoom 
useful for investigation of fine details of the iterals. 
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